Photon storage in A-type optically dense atomic media. I. Cavity model 
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In a recent paper [Gorshkov et at, Phys. Rev. Lett. 98, 123601 (2007)], we used a universal 
physical picture to optimize and demonstrate equivalence between a wide range of techniques for 
storage and retrieval of photon wave packets in A-type atomic media in free space, including the 
adiabatic reduction of the photon group velocity, pulse-propagation control via off-resonant Raman 
techniques, and photon-echo-based techniques. In the present paper, we perform the same analysis 
for the cavity model. In particular, we show that the retrieval efficiency is equal to C/(l + C) 
independent of the retrieval technique, where C is the cooperativity parameter. We also derive the 
optimal strategy for storage and, in particular, demonstrate that at any detuning one can store, 
with the optimal efficiency of C/(l + C), any smooth input mode satisfying TC7 3> 1 and a certain 
class of resonant input modes satisfying TC7 ~ 1, where T is the duration of the input mode and 2"/ 
is the transition linewidth. In the two subsequent papers of the series, we present the full analysis 
of the free-space model and discuss the effects of inhomogeneous broadening on photon storage. 

PACS numbers: 42.50.Gy, 03.67.-a, 32.80.Qk, 42.50.Fx 



I. INTRODUCTION 



The faithful storage of a traveling light pulse in an 
atomic memory and the subsequent retrieval of the state 
are currently being pursued in a number of laboratories 
around the world 0, 0, |, 1, i, i, 0, S, & M, HO, M 
US UH [HI ■ A strong motivation for this research comes 
from the field of quantum communication, where quan- 
tum information is easily transmitted by photons, but 
the photonic states need to be stored locally to process 
the information. Such applications as well as other ideas 
from quantum-information science have led to a strong 
interest in techniques to facilitate a controlled interaction 
between atoms and single photons [H, [ljj • A conceptu- 
ally simple realization of a matter-light quantum inter- 
face consists of a single atom absorbing a single photon. 
However, due to the very weak coupling of a single atom 
to light, this approach is extremely challenging and re- 
quires the use of very high-finesse cavities to effectively 
increase the coupling []J, |2fl. To circumvent the problem 
of weak coupling, it has recently been realized that one 
can use an optically thick ensemble of atoms, and sev- 
eral different proposals have been made for how external 
classical control fields can be used to controllably map 
photon states onto collective atomic states [!, 0, [H, |2C| . 
The goal in all of these approaches is to map an incom- 
ing signal pulse into a long-lived atomic coherence (re- 
ferred to as a spin wave) , so that it can be later retrieved 
"on demand" with the highest possible efficiency. Re- 
markable experimental progress has already been made 
toward the implementation of these protocols in atomic 
gases [1, 0, H, H, 0, lES and in impurities embedded in a 
solid state material [lOl [Til H2I [T3L [l4| . A central question 
that emerges from these advances is which approach rep- 
resents the best possible strategy for given experimental 
parameters and for desired memory characteristics, and 



how the control fields or possibly the shape of the input 
photon wave packet can be chosen to achieve the max- 
imum efficiency. In a recent paper (2lj j . we presented a 
novel physical picture that unifies a wide range of differ- 
ent approaches to photon storage in A-type atomic media 
and yields the optimal control strategy. This picture is 
based on two key observations. First, we showed that 
the retrieval efficiency of any given stored spin wave de- 
pends only on the optical depth of the medium and not 
on the properties of the control pulse. Physically, this 
follows from the fact that the branching ratio between 
collectively enhanced emission into desired modes and 
spontaneous decay depends only on the optical depth. 
The second observation is that the optimal storage pro- 
cess is the time reverse of retrieval (see also [2(J US H3 ) ■ 
This universal picture implies that the maximum effi- 
ciency for the combined process of storage followed by 
retrieval is the same for all approaches considered and 
depends only on the optical depth [24(. The optimum 
can be attained by adjusting the control or the shape of 
the photon wave packet. In the present paper and in the 
two papers that follow, Refs. [25|, l26| , which we will re- 
fer to henceforth as paper II and paper III, respectively, 
we present all the details behind this universal picture 
and the optimal control shaping that it implies, as well 
as consider several extensions of this analysis beyond the 
results of Ref. [21]. In particular, in the present paper 
we discuss the cavity model to be compared in paper II 
to the free-space model. In paper II, the full analysis of 
the free-space model is presented, and, in addition, the 
effects of spin-wave decay and of nondegeneracy of the 
two lower levels of the A system are discussed. Finally, 
in paper III, we generalize our treatment to two different 
regimes of inhomogeneous broadening: with and with- 
out redistribution between frequency classes during the 
storage time. 
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FIG. 1: (Color online) A-type medium coupled to a classical 
field (solid) with Rabi frequency Q(t) and a quantum field 
(dashed). Due to collective enhancement [27j ] . the quantum 
field couples to the spin-wave excitations with an effective 
coupling constant gy/~N, where N is the number of atoms in 
the medium. 



A generic model for a quantum memory uses the A- 
type level configuration shown in Fig. [H in which a weak 
(quantum) signal field (the dashed line) is detuned by 
a frequency A from the \g) — |e) transition, whose op- 
tical coherence decays with rate 7 > 7e/2, where 7 e 
is the spontaneous emission rate from state |e). A co- 
propagating (classical) control beam (the solid line) with 
the same detuning A from the \s) — |e) transition and 
time-dependent Rabi frequency envelope Q is used to 
coherently manipulate the signal propagation and map 
the photonic state onto the atoms, and vice versa. In 
the present paper and in papers II and III, we discuss 
several different approaches to photon storage, includ- 
ing far-off-resonant Raman, electromagnetically induced 
transparency (EIT), and photon-echo techniques. If we 
neglect the decay of the |s) — \g) coherence, i.e., the decay 
of the spin wave, the only sources of loss in all of these 
approaches are the decay 7 of the optical polarization on 
the \g) — |e) transition during both storage and retrieval, 
and the leakage of the pulse through the medium during 
storage. To achieve the maximum storage efficiency, one 
has to minimize these two types of loss, and, in fact, as 
we will show in the present paper and in papers II and 
III, one has to make a compromise between them. 

Higher optical depth increases the coherent coupling 
between the quantum signal field and the atoms and, 
thus, allows for higher photon storage efficiencies. It 
has therefore been su gges ted to put a cavity around the 
atomic ensemble (lTl |28|. l29j . which will effectively in- 
crease the free-space optical depth d by a factor equal to 
the number of passes a photon makes in the cavity be- 
fore leaking out. We will denote this increased effective 
optical depth by the cooperativity parameter C. High- 
efficiency retrieval of a photon from an ensemble enclosed 
in a cavity has been recently demonstrated @. In ad- 
dition to being a promising experimental setup in itself 
[30| , the slightly simpler theoretical treatment of the cav- 
ity model offers a very useful tool for understanding the 
more complicated free-space model. Thus, in the present 



paper, we will treat photon storage in the cavity model, 
to be compared in paper II to the free-space model. 

We will now review the three photon storage pro- 
tocols (Raman, electromanetically induced tranparency, 
and photon echo) that arc often discussed in the litera- 
ture on photon storage and that we treat as special cases 
of our general formalism. The remainder of this section 
is intended as an introduction to both the present pa- 
per and paper II and will thus make use of the figures of 
merit of both the cavity model (the cooperativity param- 
eter C) and the free-space model (the optical depth d). 
It will be implied in the following discussion that all the 
formulas containing C (d) refer to the cavity (free-space) 
model. 

One possible strategy for light storage uses the Ra- 
man configuration, where the fields have a large detuning 
(we will show that the appropriate limit is |A| 3> ^d or 
|A| ^> 7C rather than |A| ^ 7, as one might naively 
assume by analogy with the single-atom case) and the 
photons are absorbed into the sta ble g round state \s) by 
stimulated Raman transitions (l8l. l3ll |32|. With far-off- 
resonant interactions, the excited state |e) can be adia- 
batically eliminated to give simplified and solvable equa- 
tions [H, |34|. Based on these simplified equations, the 
Raman scheme for storage of quantum states of light in 
atomic ensembles was proposed in Ref. [lU and, simul- 
taneously with the present work, has been analyzed in 
detail and optimized under the constraint of limited con- 
trol power [3l| . We show in the present paper and in 
paper II that, in the limit of large cooperativity param- 
eter C or large optical depth d, one can ignore the d ecay 
7 of the optical \g) — |e) coherence, as is done in Ref. [3l| . 

An alternative storage strategy is based on electromag- 
netically induced transparency [H, H3, [H[ , where reso- 
nant control fields (|A| <c 'yd or |A| -C 7C) are used to 
open a spectral transparency window for the quantum 
field. In this approach, the quantum field travels at a re- 
duced group velocity, which is then adiabatically reduced 
to zero. Similarly to the Raman case, the excited state 
can also be eliminated on resonance, provided the con- 
trol field is sufficiently weak. This again simplifies the 
equations to analytically solvable ones (3f|. 

We will treat both the far-off-resonant Raman scheme 
and the resonant EIT scheme as special cases of a more 
general "adiabatic" limit, i.e., the limit in which the ex- 
cited state can be adiabatically eliminated. We will show 
that, for the purposes of optimal photon storage, the 
condition of validity of the adiabatic elimination is al- 
most independent of the single-photon detuning A (in 
particular, it is similar in the Raman and resonant cases) 
and is given by Tdj > 1 or TCj 3> 1, where T is the 
duration of the incoming pulse. We will show that, pro- 
vided a (smooth) incoming photon wave packet is long 
enough that this condition is satisfied (T ^ 1/(^7) or 
T 1/(6*7)), it can be stored with the maximum pos- 
sible efficiency, which depends only on the optical depth 
d or the cooperativity parameter C and not on the de- 
tuning A or the shape of the wave packet. In the case 
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of the cavity model discussed in the present paper, this 
maximum efficiency is simply given by C/(l + C). 

Finally, in the photon-echo-based approach to storage, 
a resonant photon is first allowed to get absorbed by the 
ensemble with the control field off. While the phrase 
"photon echo" often refers to a wide class of experiments, 
we shall here consider a special case where one applies 
then a short resonant 7r pulse, which maps excitations 
from the unstable excited state |e) into the stable ground 
state | a). Because this approach uses very short control 
pulses and because, as we will show, it is most efficient in 
storing short input pulses (T ~ 1/(^7) or T ~ 1/(C^)), 
we will refer to this approach as "fast" st orag e. This 
technique was originally suggested in Rcf. [20] for the 
case of Doppler-broadened atoms and has since been ex- 
tensively studied both theoretically [l5l I22I [H, [26| and 
experimentally [TUHH- In Ref. [2(|, it was noted that, if 
the photons are retrieved by using a control laser pulse 
traveling in the backward direction compared to storage, 
the Doppler broadening is reversed, and the dephasing 
occurring during storage is therefore also reversed, re- 
sulting in an "echo," which may result in high efficien- 
cies. In Refs. [HI HH, it was proposed to use controlled 
reversible inhomogeneous broadening (CRIB), that is, to 
artificially add an inhomogeneous broadening to an orig- 
inally homogeneously broadened line and then to reverse 
this broadening to achieve an echo signal. In the present 
paper and in paper II, we consider a different limit of this 
proposal, where there is no inhomogeneous broadening of 
the optical transition, and storage is simply achieved by 
applying a fast ir pulse at the right time. Retrieval, which 
is accomplished with a second it pulse, results [36jj in a di- 
rectional output (as opposed to the loss due to the decay 
rate 7) exactly as in the adiabatic limit. We will show 
that in the limit of large d or C this procedure leads 
to an ideal storage and retrieval of the photonic state, 
while at every finite value of d or C there exists an opti- 
mal input photon mode that can be stored with efficiency 
equal to the maximum adiabatic storage efficiency (given 
by C/(l + C) in the case of the cavity model). For com- 
parison, in paper III, we will discuss how this approach 
measures up to the CRIB approach and show that adding 
and reversing inhomogeneous broadening as proposed in 
Refs. [HI, H3| may lead to an improvement in the storage 
efficiency, although the improvement is rather limited. 

The optimization of storage in all of these schemes con- 
sists of finding the optimal balance between two sources 
of error: leakage of the input pulse through the ensemble 
and scattering of the input photons into Att due to sponta- 
neous emission. In the EIT approach, a stronger control 
field is desirable, since it produces more robust interfer- 
ence and a wider transparency window, thus minimizing 
spontaneous emission losses. On the other hand, higher 
control power means larger group velocity and hence the 
inability to localize the input pulse inside the medium. 
The optimization in this case finds the optimal power and 
shape for the control field, given the duration and shape 
of the input pulse. In contrast, in the Raman scheme, a 



high value of fl is required to have a sufficient coupling 
of the input photon to the spin wave [3lj . On the other 
hand, large Q will increase the decay rate due to spon- 
taneous emission, which is given by the optical pumping 
rate j£l 2 / A 2 . The optimization with respect to the shape 
and power of f2 for a given input mode again balances 
between these two sources of error. Finally, in the fast 
storage scheme the control fields are fixed to be perfect 
7r pulses, but one can optimize with respect to the dura- 
tion T and the shape of the input mode. The input mode 
should be made as short as possible in order to avoid the 
loss due to optical polarization decay exp(— 7T). How- 
ever, a mode that is too short will be too wide in fre- 
quency space and will not be absorbed by the ensemble 
(i.e., it will leak through). The optimization with respect 
to the duration and shape of the input mode finds the 
optimal balance between these two sources of error. 

In all the photon storage techniques considered, ideal 
performance (i.e., unit efficiency) can be achieved in the 
limit of infinite optical depth d or infinite cooperativity 
parameter C . For example, in the EIT regime in free 
space, the width of the spectral transparency window is 
Aweit = VgVd/L, where L is the length of the ensemble 
and u g oc \fl\ 2 /d is the EIT group velocity [37j]- Thus, 
for a given T and a given large value of d, one can first 
make Q, and hence v g , small enough for the pulse to 
fit inside the medium. Then the enhancement of Albeit 
by an extra factor of \fd will ensure, if d is sufficiently 
large, that the transparency window is still wide enough 
to induce negligible spontaneous emission. In the Raman 
regime, to avoid spontaneous emission decay via the opti- 
cal pumping rate 7i7 2 /A 2 , one should make Q, sufficiently 
small. If the optical depth or the cooperativity parame- 
ter is large enough, the coupling of the input mode to the 
atoms will then still be sufficient to avoid leakage even at 
this small value of tt. Finally, in fast storage, the pulse 
that is short enough (T7 <C I) to avoid optical polariza- 
tion decay can still be absorbed in a free-space medium 
provided d is large enough (Tdj ^> I, as we will show 
in paper II). In the cavity model discussed in the present 
paper, due to the availability of only one spin-wave mode 
(the one that couples to the cavity mode), high-efficiency 
fast storage is harder to achieve than in free space: only 
pulses of a particular shape and duration (T ~ I/(C7)) 
give high fast storage efficiencies. 

Although ideal performance can be achieved at infi- 
nite optical depth, in practice, optical depth is always 
limited by experimental imperfections such as a limited 
number of atoms in a trap (e.g., Ref. [?}, where the opti- 
cal depth is roughly 8), competing four- wave mixing pro- 
cesses (e.g., Ref. p, where the optical depth is roughly 
4), inhomogeneous broadening of impurity levels in solid 
state samples 15[ , or other types of experimental imper- 
fections. Therefore, the optimization of storage protocols 
at finite optical depth is essential. 

Before proceeding with our analysis, we would like to 
specially note the recent work of Dantan et at, which 
also considers and illuminates some of the issues we dis- 
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cuss in the present paper [38|, 133. l4Ci |41| and in paper 
II [ilj]. In particular, in Refs. |38l. l39l. l40j|. focusing on 
broadband squeezed states as the input, the authors con- 
sider adiabatic storage in a cavity, derive an efficiency ex- 
pression equivalent to ours, and recognize the interesting 
similarity between Raman and resonant regimes, both of 
which feature reduced sensitivity to spontaneous emis- 
sion. We show in the present paper how, for the case 
of a single incoming spatiotemporal field mode, proper 
control field shaping can be used to achieve the same op- 
timal efficiency independent of detuning. This effectively 
makes the Raman, the resonant, and the intermediate 
regimes all exhibit equally reduced sensitivity to sponta- 
neous emission. In Ref. [41] . as in paper II, the authors 
compare adiabatic storage in a cavity to adiabatic stor- 
age in free space and recognize important similarities and 
differences. In particular, it is recognized that the error 
in the cavity case scales as the inverse of the atomic den- 
sity (as we also find in the present paper), while in free 
space it may scale as the inverse of the square root of the 
density. We show in paper II how proper optimization 
can be used to make the error in the free-space case also 
scale as the inverse of atomic density. We also bridge in 
free space the gap between the EIT and Raman cases, 
showing how proper control field shaping can be used to 
achieve the same optimal efficiency independent of de- 
tuning. 

We would also like to note that the connection be- 
tween optimal photon storage and time reversal that we 
present was first made in the context of photon-echo- 
based techniques. In particular, it was first shown in 
Ref. [2(| and then discussed in detail in Refs. [IS HH that 
under certain conditions, such as high optical depth and 
sufficiently slow optical polarization decay rate, photon- 
echo techniques can result in ideal storage and retrieval, 
and that the retrieved photon field is then the time re- 
verse of the original input field. We generalize this result 
in Ref. [2l|, in the present paper, and in papers II and III 
by demonstrating that the ideas of time reversal can be 
used to optimize photon storage even when the dynamics 
of the system are not completely reversible and when the 
ideal unit efficiency cannot be achieved. This is the case 
for finite cooperativity parameter C in the cavity model 
and for finite optical depth d in the free-space model. 
We also generalize the time-reversal-based optimization 
of photon storage from photon-echo-based techniques to 
any storage technique including, but not limited to, EIT 
and Raman techniques in homogeneously (present paper, 
Ref. [21| , and paper II) and inhomogeneously (paper III) 
broadened A-type media. 

We would also like to point out that mathematically 
some of the optimization problems we are solving in the 
present work (including Ref. [2l|, the present paper, and 
papers II and III) fall into a rich and well-developed field 
of mathematics called optimal control theory [42l . |43| . |44| . 
In particular, we are interested in shaping the control 
pulse to maximize the storage efficiency (or efficiency 
of storage followed by retrieval) for a given input pho- 



ton mode in the presence of optical polarization decay. 
Since in the equations of motion the control multiplies 
a dependent variable (optical polarization) , this problem 
is a nonlinear optimal control problem [42] , A general 
solution to all nonlinear optimal control problems does 
not exist, so that such problems have to be treated on 
a case by case basis. We believe that the methods we 
suggest, such as the iterative time-reversal method in- 
troduced in Sees. IV and V of paper II, may be use- 
ful in solving optimal control problems in other open 
(as well as closed) quantum systems. Similar iterative 
methods are a standard tool in applied optimal control 
US EE EE El] and have been used for a variety of ap- 
plications, including laser control of chemical reactions 
HE El , design of NMR pulse sequences [49| , loading of 
Bose-Einstein condensates into an optical lattice 15 CM . and 
atom transport in time-dependent superlattices [5l|. In 
fact, an optimization procedure that is based on gradi- 
ent ascent 1491 . and that is very similar to that of Refs. 
[H El, EEHfl E3, EE El, is directly applicable to 
our problem of finding the optimal control pulse, as we 
will discuss elsewhere [52] . However, in the present paper 
and in papers II and III, we use time-reversal iterations 
for optimal control in a way different from the methods of 
Refs. 0, El, EE EE ES EE EI HE HH , as we will discuss 
in Sec. V of paper II. In particular, we will show that, in 
addition to being a convenient computational tool, our it- 
erative optimization algorithm is, in fact, experimentally 
realizable (53|. 

The remainder of this paper is organized as follows. 
Section [II] applies to both the cavity and the free-space 
models and discusses our figure of merit for the perfor- 
mance of the photon storage. The rest of the paper dis- 
cusses storage and retrieval of photons using homoge- 
neously broadened atomic ensembles enclosed in a cav- 
ity. In Sec. IIIH we introduce the model. In Sec. HVl 
without fully solving the equations analytically, we show 
that both the retrieval efficiency and the optimal stor- 
age efficiency are equal to C/(l + C) (where C is the 
cooperativity parameter), and derive the optimal stor- 
age strategy. In Sees. [V] and IVI1 we solve the equations 
analytically in the adiabatic and fast limits, respectively, 
and demonstrate that the optimal storage efficiency can 
be achieved for any smooth input mode at any detuning 
satisfying TC7 3> 1 and a certain class of resonant input 
modes satisfying TC7 ~ 1, where T is the duration of the 
input mode. In Sec. IVII1 we summarize the discussion of 
the cavity model. Finally, in the Appendixes, we present 
some details omitted in the main text. 



II. FIGURE OF MERIT 

When comparing different storage and retrieval ap- 
proaches, it is essential to have a figure of merit charac- 
terizing the performance of the memory. The discussion 
in this section of the appropriate figure of merit applies 
both to the cavity model discussed in this paper and 
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to the free-space models discussed in papers II and III. 
Throughout this work we shall assume that we initially 
have a single incoming photon in a known spatiotem- 
poral mode denoted by £- m (t) (or, for the case of com- 
puting retrieval efficiency alone, a single excitation in a 
known atomic spin- wave mode). We define the efficiency 
rj of all the mappings we consider (storage alone, retrieval 
alone, or storage followed by retrieval) as the probability 
to find the excitation in the output mode (photonic or 
atomic, as appropriate) after the interaction. Depending 
on the application one has in mind, this single-photon 
efficiency may or may not be the right quantity to con- 
sider, but provided that we are interested in a situation 
where we are mapping a single input mode into a single 
output mode, any other quantities may be derived from 
the single-photon efficiency rj. 

For all the interactions we consider, the full evolution 
results in a passive (beam-splitter-like) transformation 

bj=^2u jk a k , (I) 

k 

where aj and b k denote the annihilation operators for all 
the input and output modes, respectively (all photonic, 
spin- wave, and Langevin noise operators), with commu- 
tation relations [dj-,dt] = Sj >k and [Sj,St] = Sjk. Here 
the matrix U has to be unitary to preserve the com- 
mutation relations. The mapping from a certain input 
mode So to an output mode bo with efficiency rj is there- 
fore described by 60 = y^cto + \/l — rjc, where c satisfies 
[c, C*} = 1 and represents some linear combination of all 
other input modes orthogonal to d . If all input modes 
other than do are in the vacuum state, the parameter rj 
completely characterizes the mapping. If, for instance, 
the mode we are storing is in an entangled state with 
some other system (|0)a o |x) + |l)a |y))/2, where |0)a o and 
|l)a are the zero- and one-photon Fock states of the in- 
put mode, and \x) and \y) are two orthonormal states of 
the other system, the fidelity of the entangled state after 
the mapping is easily found to be F = (1 + rj)/2. Simi- 
larly, Refs. [38|, [33, Ho, [IjJ characterize the performance 
in terms of squeezing preservation parameter rj squcczc . If 
the input state is a squeezed vacuum state in a given 
mode do, the squeezing preservation parameter can be 
shown to be equivalent to single-photon efficiency, i.e., 

^squeeze 

= rj. We will show below in the description of our 
model why in most experimental situations it is indeed 
reasonable to assume that the incoming noise (which is 
included in c) is vacuum noise. 



III. MODEL 

The details of the model and the derivation of the equa- 
tions of motion are provided in Appendix [X] In this sec- 
tion, we only give a brief introduction to the model and 
present the equations of motion without derivation. 



We consider a medium of N A-type atoms with two 
metastable lower states, as shown in Fig. [TJ interacting 
with two single-mode fields. We neglect reabsorption of 
spontaneously emitted photons and treat the problem in 
a one-dimensional approximation. The \g) — |e) transi- 
tion of frequency uj eg of each of the atoms is coupled to 
a quantized traveling- wave cavity radiation mode (e.g., a 
mode of a ring cavity with one of the mirrors partially 
transmitting) with frequency u>i = uj eg — A described by 
a slowly varying annihilation operator £{t). The cavity 
decay rate is 2k and the corresponding input-output re- 
lation is [H3 | 

£ out (t) = V2k£(t)-£ in {t). (2) 

In addition, the transitions \s) — |e) of frequency lu £S are 
driven by a single-mode copropagating classical plane- 
wave control field with frequency L02 = u) es — A (i.e., 
at two-photon resonance u)\ — UJ2 = to S g, where Hu) sg is 
the energy difference between the two metastable states) 
described by a slowly varying Rabi frequency envelope 

In the dipole and rotating-wave approximations, as- 
suming that almost all atoms are in the ground state at 
all times, and defining the polarization annihilation op- 
erator P(t) = a ge (t)/\/~N and the spin- wave annihilation 
operator S(t) — a gs (t)/ yN (where a^ v are slowly vary- 
ing collective atomic operators defined in Appendix |A"|) . 
to first order in £ , the Heisenberg equations of motion 
are 

£ = -K£ + ig^P + V2K~£ in , (3) 
P = -{-f + iA)P + igVN£ + tnS+^F P , (4) 
i = -%S + iQ*P+y/2T s F s , (5) 

where we have introduced the polarization decay rate 
7, the spin- wave decay rate j s , and the corresponding 
Langevin noise operators Fp and F$- The coupling con- 
stant g (assumed to be real for simplicity) between the 
atoms and the quantized field mode is collectively en- 
hanced H?! by a factor of \/N to g\f~N ■ 

As described in Appendix IA1 under reasonable exper- 
imental conditions, the incoming noise described by Fp 
and Fs is vacuum, i.e., all normally ordered noise cor- 
relations are zero. This is precisely the reason why, as 
noted in Sec. UH efficiency is the only number we need in 
order to fully characterize the mapping. 

We assume that all atoms are initially pumped into the 
ground state, i.e., no P or S excitations are present in 
the atoms. We also assume that the only input field ex- 
citations initially present are in the quantum field mode 
with an envelope shape ho(t) nonzero on [0,T]. The goal 
is to store the state of this mode in S and, starting at 
a time T r > T, retrieve it back into a field mode. Since 
we are interested only in computing efficiencies (defined 
below) and since the incoming noise is vacuum, we can 
ignore the noise operators in Eqs. ([2])-© and treat these 
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equations as complex number equations. During stor- 
age, the initial conditions are P(0) = 0, 5(0) = 0, and 
the input mode is £- ln {t) = ho(t) (normalized according 
to J dt\£ in (t)\ 2 = 1). We have here dropped the carets 
on the operators to denote their complex number repre- 
sentations. The storage efficiency is then 



Va = 



(number of stored excitations) 
(number of incoming photons) 



\S(T)\ 2 . (6) 



Similarly, during retrieval, the initial and boundary con- 
ditions are P(T r ) = 0, S(T T ) = S(T), and £ in (t) = 0. 
£out{t) then represents the shape of the quantum mode 
into which we retrieve, and the total efficiency of storage 
followed by retrieval is given by 



Vtot 



(number of retrieved photons) 
(number of incoming photons) 



dt\£ out (t)\< 



(7) 

If we instead take S(T T ) = 1, we obtain the retrieval 
efficiency: 



(number of retrieved photons) 
(number of stored excitations) J T 



dt\£ out (t)\ 2 



(8) 

From now on we will neglect the slow decay of the spin 
wave (i.e., set 7 S = 0) but, as briefly discussed below 
at the ends of Sees. IV Al and IV Bl spin- wave decay is 
not hard to include. Nonzero 7 S will simply introduce 
an exponential decay without making the solution or the 
optimal control shaping harder. 

To get the closest analogy to the free-space regime, 
we assume we are always in the "bad cavity" limit 
(k gyN), in which £ in Eq. ((3]) can be adiabatically 
eliminated to give 



£o 



£ m + iy/2jCP, 



(9) 



P = -(7(1 + C) +iA)P + inS + i v / 2^fC£ in ,{10) 

s = m*p, (11) 

where C — g 2 N / (k'j) is the cooperativity parameter. 
To relate to the free space situation discussed in pa- 
per II, we can write the cooperativity parameter as 
C = 2d[(l/(2/s))/(L/c)] ) where d = g 2 NL/(cj) is the 
definition of optical depth used in the free-space model 
of paper II and where the factor in the square brackets 
(proportional to cavity finesse) is equal to the number 
of passes a photon would make through an empty cavity 
before leaking out (i.e., the photon lifetime in the cavity 
divided by the time a single pass takes). Thus, up to a 
factor of order unity, the cooperativity parameter C rep- 
resents the effective optical depth of the medium in the 
cavity, so that the efficiency dependence on C in the cav- 
ity should be compared to the efficiency dependence on 
d in free space. We note that, although Eqs. (j9])- (fTTj) de- 
scribe our case of quantized light coupled to the \g) — |e) 
transition, they will also precisely be the equations de- 
scribing the propagation of a classical probe pulse. To 



see this one can simply take the expectation values of 
Eqs. ©-© and use the fact that classical probe pulses 
are described by coherent states. 

It is convenient to reduce Eqs. (fTU|) and (fTTj) to a single 
equation 



n* 



S 



- (7(1 + C) + iA)S+ \n\ 2 S = -fi* V27C£in- 

This second-order differential equation cannot, in gen- 
eral, be fully solved analytically. However, in the next 
section we derive a number of exact results about the 
optimal efficiency anyway. 



IV. OPTIMAL STRATEGY FOR STORAGE 
AND RETRIEVAL 

In this section, we derive several important results re- 
garding the optimal strategy for maximizing the storage 
efficiency, the retrieval efficiency, and the combined (stor- 
age followed by retrieval) efficiency without making any 
more approximations. 

It is convenient to first consider retrieval. Although we 
cannot, in general, analytically solve for the output field 
£out(t), we will now show that the retrieval efficiency is 
always C/(l + C) independent of the detuning A and the 
control shape £l(t) provided that no excitations are left 
in the atoms at t = 00, i.e., P(oo) = and S(oo) = 0. 
From Eqs. (JTQJ) and (fTTj) and using £i n (t) = 0, we find 



dt 



(\P\ 2 + \S\ 2 )=-2j(l + C)\Pf 



(13) 



Using this and Eqs. © and ([5]), the retrieval efficiency 
becomes 



n 



c 



(\S(T T )\ 2 + \P(T X )\* -\S(oc)\ 2 -\P(<x>)f 



l + C 

(14) 

which reduces to C/(l + C) for S{T X ) = 1, P(T r ) = 
P(oo) = 5(oo) = 0. The value of the retrieval error 
(1 — rfi = 1/(1 + C)) and its independence from A and 
il follow directly from the branching ratio between the 
decay rates in Eq. (fT0|) (or equivalently in Eq. (fl3|) ). The 
decay rate for P into undesired modes is 7, while the 
decay rate for P into the desired mode £ ou t is 7C. The 
retrieval efficiency, which is the ratio between the desired 
decay rate and the total decay rate, is, therefore, equal 
to C/(l + C) independent of the control field. 

We have thus shown that, provided our control pulse is 
sufficiently long and/or powerful to leave no excitations 
in the atoms (we will refer to this as complete retrieval) , 
the retrieval efficiency is independent of A and f2(t) and 
is always equal to CJ (l + C). Therefore, any control field 
is optimal for retrieval provided it pumps all excitations 
out of the system. Using this knowledge of the retrieval 
efficiency, in the remainder of this section we will use 
a time-reversal argument to deduce the optimal storage 
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strategy and the optimal storage efficiency. Here we will 
only give the essence of and the intuition behind the time- 
reversal argument, and leave the derivation to Sees. IV 
and V of paper II. In the remainder of the paper, we will 
independently confirm the validity of this argument in 
the adiabatic and fast limits. 

Applied to the present situation, the essence of the 
time-reversal argument is as follows. Suppose one fixes 
the cooperativity parameter C and the detuning A and 
considers complete retrieval from the spin wave with a 
given control field Cl(t) into an output mode £ ou t(t) of 
duration T ou t. According to the time-reversal argument, 
the efficiency for storing the time reverse of the output 
field (£ in (t) = £ * ut (T 0Ut - *)) with n*(T ont - t), the time 
reverse of the retrieval control field, into the spin wave 
is equal to the retrieval efficiency [55J. Although this 
claim is not trivial to prove (see paper II), it is rather 
intuitive: since the retrieval procedure can be regarded as 
a generalized beam-splitter-like transformation (Sec. HT| . 
the equality of the two efficiencies is simply the statement 
that the probability of going from a given input port of 
the beam splitter to a given output port is equal to the 
probability of going backward from that output port to 
the original input port. 

Therefore, the time-reversal argument shows that the 
maximum efficiencies for storage and storage followed by 
retrieval are C/(l + C) (i.e., the retrieval efficiency) and 
C 2 /(l + C) 2 (i.e., its square), respectively. Moreover, 
it says that these maximum efficiencies are obtained if 
the input field £ m (t) and the storage control field f2(i) 
are such that fl*(T — t), i.e., the time reverse of f2(t), 
retrieves the spin-wave excitation into the output mode 
£out(i) = £in(T — t), i.e., the time reverse of £i n (t). I n 
order to say for which input fields the optimal storage 
control (l(t) can be found (or, equivalently, into which 
output fields a spin- wave excitation can be retrieved), 
we need to consider the limits, in which Eq. (|12p can be 
fully solved analytically. These limits, adiabatic and fast, 
will be discussed in the following sections. 



V. ADIABATIC RETRIEVAL AND STORAGE 

A. Adiabatic retrieval 

In the previous section, we have found, based on time 
reversal, the maximum storage efficiency and the scenario 
under which it can be achieved. Since the optimal stor- 
age into a given input mode requires the ability to carry 
out optimal retrieval into the time reverse of this mode, 
we will, in the following sections fSecs. FVl and lVI[) . solve 
Eq. (|12[) analytically in two important limits to find out 
which modes we can retrieve into and store optimally. 
The first such limit, which we will discuss in this section 
(Sec. |V|, corresponds to smooth control and input fields, 
such that the term in the square brackets in Eq. (fT2")l can 
be dropped. This "adiabatic" limit corresponds to an 
adiabatic elimination of P in Eq. (JTO|) . The precise condi- 



tions for this approximation will be discussed in Sec. IV Cl 
In this section, we discuss the retrieval process. 

It is instructive to recognize that in the adiabatic ap- 
proximation (i.e., with P in Eq. (fTO)) replaced with 0), if 
one uses rescaled variables £i n (t)/il(t), £ out (t)/il(t), and 
P(t)/Q(t) and makes a change of variables t — > h(T T ,t), 
where 



h(t,t') 



\n(t")\ 2 dt'' 



(15) 



then Eqs. (j!?j)- (fTTj) become independent of f2 and can be 
solved in this ^-independent form, so that for any given 
fi the solution in the original variables would follow by 
simple rescaling. However, since the equations are suffi- 
ciently simple and in order to avoid confusion introduced 
by additional notation, we will solve Eqs. (f9|)- (fTTj) directly 
without making the change of variables. 

To compute the output field during adiabatic retrieval, 
we assume for simplicity that retrieval begins at time t = 
rather than at time t = T T and adiabatically eliminate 
P in Eqs. (JTUJ) and (JTTJ) (i.e., replace P in Eq. (TO]) with 
zero) to obtain a first-order linear ordinary differential 
equation for 5. Then, using 5(0) = 1 and £- m (t) = 0, we 
solve this equation to find 



£ont(t) = -V2^C 



Q(t) 



7 (l + C*) + iA 



e t (i+c)+;a 



h(0,t) 



(16) 



The ^-dependent phase ih(0, t)A/(j 2 (l + C f + A 2 ) in 
the last factor is the ac Stark shift, which results in a 
shift of the output field frequency away from bare two- 
photon resonance. Computing the retrieval efficiency us- 
ing Eq. (|16p. we find 



Vr 



c 



1 + C 



27(1+0) 
1 _ e ^(i+C^+a: 



-h(0,oo) 



(17) 



which is equal to C /{1 + C) provided the control pulse is 
sufficiently powerful and/or long to ensure that 



2 7 (1 + C) 
7 2 (1 + C) 2 + A 2 



/i(0,oo) > 1, 



(18) 



which is the same as the condition P(oo) = 5(oo) = 0. 
Note that adiabatic elimination did not affect the exact 
value of the efficiency and kept it independent of fl(t) and 
A by preserving the branching ratio between the desired 
and undesired state transfers. Also note that, unlike the 
general argument in the previous section, which assumed 
P(oo) = 5(oo) = 0, Eq. (fTT)) allows for the precise calcu- 
lation of the retrieval efficiency for any /i(0, oo). 

As noted in the Introduction, two important subsets 
of the adiabatic limit, the resonant limit and the Raman 
limit, are often discussed in the literature. Although, 
as we show in this work, the basic physics based on the 
branching ratio and time-reversal arguments is shared by 
both of these approaches to quantum memory, a more 
detailed discussion of the physics behind them involves 
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significant differences. In fact, prior to this work, the fact 
that the two approaches are in a sense equivalent was not 
recognized to our knowledge: only interesting similarities 
were pointed out (38l . l39l l4l| . As an example of an im- 
portant difference, the resonant and Raman limits give 
different dependences on C of the duration of the output 
pulse in Eq. IT6|) : 



j 2 C 2 + A 2 



(19) 



where we assumed C > 1. In the resonant limit (7C 
|A|), r ou t ~ jC/\Q\ 2 , while in the Raman limit (7C < 
|A|), T out ~ A 2 /(7C|rj| 2 ). It is worth emphasizing that 
the Raman limit condition is 7C -C |A| and not 7 -C |A|, 
as one may naively think by analogy with the single-atom 
case. 

It follows from the concept of time reversal that the 
modes that can be stored optimally are the time reverses 
of the modes onto which a spin wave can be retrieved. 
We will now show that, in the adiabatic limit, at any 
given A and C, we can shape f2(t) to retrieve onto any 
normalized mode e(t). Integrating the norm squared of 
Eq. HID with £ out (t) = y/C/(l + C)e(t), we get 



dt'\e(t')\ 2 = 1-e 



(20) 



Solving this equation for /i(0, t) and then taking the 
square root of its derivative with respect to t, we find 
|fi(t)|. Knowing h(0,t), the phase of Q(t) can be deter- 
mined from Eq. (|16|) . Putting the magnitude and the 
phase together, we have 



7(1 + C) + iA 



it) 



Ah(Q,t) 



V27(l + C) ^/ t °° \ e {t')\ 2 dt> 

(21) 

where h(0,t) should be determined from Eq. (|20| . For 
any e(t), this expression gives the control Q(i) that re- 
trieves the spin wave into that mode. The phase of Q(t), 
up to an unimportant constant phase, is given by the 
phase of the desired output mode plus compensation for 
the Stark shift (the last factor). It is also worth noting 
that, up to a minus sign and a factor equal to the first 
fraction in Eq. (|2"Tj) , Q(t) is simply equal to e(t)/S(t). 

We note that, if one wants to shape the retrieval into 
a mode e(t) that drops to zero at some time T out suffi- 
ciently rapidly, |f2(i)| in Eq. (|2"Tj) will go to 00 at t — T out . 
The infinite part can, however, be truncated without sig- 
nificantly affecting the efficiency or the precision of e(t) 
generation. One can confirm that the loss in efficiency is 
small by inserting into the adiabatic solution in Eq. (JT7J) 
a value of /i(0, 00) that is finite but large enough to satisfy 
Eq. (fl"8"|) . One can similarly confirm that the generation 
of e(t) can be precise with truncated control fields by us- 
ing Eq. (|16p . However, to be completely certain that the 
truncation is harmless, one has to solve Eqs. (f9"j)- (fTTj) nu- 
merically without making the adiabatic approximation. 



We will do this in Sec. IV Cl for the case of storage, where 
the same truncation issue is present. 

We briefly mention that the spin- wave decay rate %, 
which we have ignored so far, simply introduces a decay 
described by exp(— j a t) into Eq. (fit)]) and, unless we re- 
trieve much faster than 1/7 S , makes retrieval efficiency 
control dependent. With nonzero 7 S , we can still shape 
retrieval to go into any mode: we shape the control us- 
ing Eq. (|2"Tj) as if there were no 7 S decay except that the 
desired output mode e(t) should be replaced with the 
normalized version of e(t) exp(7 s i), i.e., 



i{t) e(t)e 7at 



dt'\e{t')\ 2 e 



(22) 



The retrieval efficiency will, however, be output-mode- 
dependent in this case: it will be multiplied (and hence 
reduced) by [ $™ dt'\e{t')\ 2 exp(2 7s i')] _1 . 



B. Adiabatic storage 

In principle, using the solution for retrieval from the 
previous section, the time-reversal argument of Sec. IV 
immediately guarantees that, provided we are in the adi- 
abatic limit (conditions to be discussed in Sec. lVCp . we 
can always shape the control field to store any input mode 
£i n (t) at any detuning A with the maximum efficiency 
(7/(1 + C). However, for completeness, and to verify 
that the optimal storage control field is indeed the time 
reverse of the control field that retrieves into £* n (T — t), 
we give in this section the solution to adiabatic storage. 

In the adiabatic approximation, we use a procedure 
very similar to that used in the retrieval solution, to find 



S(T) = 



C 



1 + C 



dtf(t)£ in (t), 



where 



0*^)^/27(1 + C) Mt.T) 

7 (l + C)+iA 



/(*) 

The storage efficiency is then 
C 



1 + C 



dtf(t)£ in (t) 



(23) 



(24) 



(25) 



We are interested in computing the control that maxi- 
mizes r] s for a given & m (t) . We find in Appendix |5] that 
the maximum storage efficiency is C / (1 + C) and that it 
can be achieved (in the adiabatic limit) for any A and 
£in(t), and that the optimal control is 



fi(t) = - 



7(l + C)-iA 



V / 27(T+C) yjf*\£ ia (t>)\*dt' 



Ah(t,T) 
-y^fl + C^+A^ 



(26) 
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where h(t, T) can be found by inserting Eq. (|26| into 
Eq. dT5j) . The phase of Cl(t), up to an unimportant con- 
stant phase, is thus given by the phase of the input mode 
plus compensation for the Stark shift (the last factor). 
As for the retrieval control discussed in the previous sec- 
tion, we note that, although |Q(t)| in Eq. ([26]) goes to oo 
at t = 0, the infinite part can be truncated without sig- 
nificantly affecting the efficiency. This can be confirmed 
analytically using Eq. (|25|) provided the adiabatic limit 
is satisfied. We will also confirm this numerically in the 
next section without making the adiabatic approxima- 
tion. 

As expected from the time- reversal argument, the op- 
timal control we derived is just the time reverse — > 
fl*(T - t)) of the control that retrieves into £* n (T - t), 
the time reverse of the input mode. We verify this in 
Appendix [B] 

Although optimal storage efficiencies are the same in 
the Raman and adiabatic limits, as in the case of re- 
trieval, rather different physical behavior can be seen in 
the two limits. It is now the dependence on C of the opti- 
mal control intensity (which can be found from Eq. (|26|) ) 
that can be used to separate resonant and Raman behav- 
ior. Assuming for simplicity C > 1, in the resonant limit 
(7(7 > |A|), |fi| - y/'yC/T, while in the Raman limit 
(7C < |A|), |0| ~ \A\/y/yVT. Since complete retrieval 
and optimal storage are just time reverses of each other, 
it is not surprising that these relations are identical to 
the ones we derived for the dependence of output pulse 
duration on C in the previous section. This opposite de- 
pendence of |f2| on C in the Raman and EIT limits is, 
in fact, the signature of a simple physical fact: while the 
coupling of the input photon to the spin wave increases 
with increasing f2 in the Raman case, it effectively de- 
creases in the EIT regime where a very large will give 
a very wide transparency window and a group velocity 
equal to the speed of light. This is why as the coopera- 
tivity parameter changes, the control has to be adjusted 
differently in the two regimes. 

As for retrieval, we briefly mention that nonzero 7 S 
simply introduces exp(— j s (T — t)) decay into Eq. 
The optimal storage control can still be found using 
Eq. (|26|) as if there were no decay, except that the input 
mode should be replaced with the normalized version of 
£ in (*)exp(-7 s (T-*)), i.e., 



dt'\£ in (t')\ 2 e 



2--2 7B (T-t') 



(27) 

However, the optimal storage efficiency will now depend 
on input pulse duration and shape: it will be multiplied 
(and hence reduced) by dt'\£ in {t')\ 2 exp(-2j a (T-tf)). 
It is important to note that with nonzero spin-wave de- 
cay the optimal storage efficiency of a particular input 
mode is no longer identical to the retrieval efficiency into 
its time reverse. This is not at variance with the time- 
reversal argument discussed in detail in paper II (which 



still applies when 7s ^ 0), since the corresponding opti- 
mal storage and retrieval control shapes are no longer the 
time reverses of each other, and, in contrast to the 7 S = 
case, the retrieval efficiency is now control dependent. 
Finally, we note that, when we consider storage followed 
by retrieval, in order to take into account the spin-wave 
decay during the storage time [T, T r ], one should just 
multiply the total efficiency by exp(— 27 s (T r — T)). 



C. Adiabaticity conditions 

We have found that, provided we are in the adiabatic 
limit, any input mode can be stored optimally. In this 
section we show that, independent of A, the sufficient 
and necessary condition for optimal adiabatic storage of 
a pulse of duration T to be consistent with the adiabatic 
approximation is TC7 1. 

To find the conditions for the adiabatic elimination of 
P in Eq. (|10[) , we do the elimination and then require its 
consistency by enforcing [56j ] 



\P\ < \(jC + iA)P\ 



(28) 



(we assume for simplicity throughout this section that 
C > 1). During retrieval, sufficient conditions for 
Eq. (HHD are 



fi| < \^C + iA\ 



n 



< \ jC + iA\, 



(29) 
(30) 



which limit, respectively, the power and the bandwidth 
of the control pulse. These are easily satisfied in practice 
by using sufficiently weak and smooth retrieval control 
pulses. 

During storage, the satisfaction of Eq. (|28|) requires, in 
addition to conditions (f2"9"| and ([50]) , the satisfaction of 



< |7<7 + «A| 



(31) 



which limits the bandwidth of the input pulse. In partic- 
ular, for a smooth input pulse of duration T, this condi- 
tion is implied by 



T<7 7 > 1. 



(32) 



Let us now show that for optimal storage, the condition 
(f3"2")l also implies conditions (|29[) and (f3"0"|) and is thus 
the only required adiabaticity condition (provided £ m is 
smooth). Application of Eq. (f2T))) reduces Eq. (|29[) to 
Eq. (f32|) . Equation (j30|) . in turn, reduces to the con- 
ditions on how fast the magnitude |fi| and the phase <f> 
(which compensates for the Stark shift) of the control can 
change: 



101 



< | 7 C + iA| 



(33) 
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where fl = |f2| exp(i(/>). Application of Eq. (|26[) shows 
that Eq. implies Eq. 

We have thus shown that TC7 3> 1 is a sufficient con- 
dition for the validity of adiabatic elimination in optimal 
storage. But, in fact, from the amplitude of the optimal 
storage control field (Eq. (j26|) ). one can see that Eq. (j29|) 
implies that TC7 3> 1 is also a necessary condition for 
the validity of adiabatic elimination in optimal storage. 
(To show that TC7 3> 1 is a necessary condition, one 
also has to use the extra condition on the adiabatic value 
of P(0) to rule out the special situation when Q(t) 
is such that Eq. (f2"5)) is satisfied but Eq. (j2"9")) is not.) 

To verify the adiabaticity condition in Eq. ([32]) and 
investigate the breakdown of adiabaticity for short input 
pulses, we consider a Gaussian-like input mode (shown 
in Fig. 3 of paper II) 



£ in (t) = ^( e -30( t /T-0.5) 2 



-7.5 



T. 



(34) 



where for computational convenience we have required 
£ in (0) = £ in (T) = and where A w 2.09 is a nor- 
malization constant. We fix the cooperativity parame- 
ter C, the detuning A, and the pulse duration T, and 
use Eq. (f2"fi|) to shape the control field. We then use 
Eqs. (|10[) and (jTTJ) without the adiabatic approximation 
to calculate numerically the actual storage efficiency that 
this control field gives, and multiply it by the control- 
independent retrieval efficiency C/(l+C), to get the total 
efficiency of storage followed by retrieval. As we decrease 
T, we expect this efficiency to fall below (C/(l + C)) 2 
once TC7 3> 1 is no longer satisfied. And indeed 
in Fig. HJa) we observe this behavior for A = and 
C = 1,10,100,1000. In Fig. Hfb), we fix C = 10 and 
show how optimal adiabatic storage breaks down at dif- 
ferent detunings A from to IOOO7. From Fig. [DJb), we 
see that, as we move from the resonant limit (C7 3> |A|) 
to the Raman limit (C7 <C |A|), we can go to slightly 
smaller values of TC7 before storage breaks down. How- 
ever, since the curves for A = IOO7 and A = IOOO7 
almost coincide, it is clear that TC7 ^ 1 is still the rel- 
evant condition no matter how large A is, which must 
be the case since the condition (|2"9")l breaks down for 
shorter T . The most likely reason why in the Raman 
limit adiabaticity is slightly easier to satisfy is because 
in the Raman limit it is only condition (|29p that reduces 
to TC7 ^> 1, while conditions (|3"0|) and (f3"Tj) reduce to 
TA » 1, which is weaker than TC7 > 1 (since A > C7 
in the Raman limit). In the resonant limit, in contrast, 
all three conditions (ESD-(EH) reduce to TC7 > 1. 

Before turning to the discussion of fast retrieval and 
storage, we note that the use of Eq. (f2"o| to calculate 
the storage control fields for Fig. [2] resulted in a control 
field f2(i) whose magnitude went to 00 at t — 0, as pre- 
dicted in the previous section. To generate Fig. [21 the 
optimal \£l(t)\ were therefore cut off for t < T/100 to 
take the value |fi(T/100)|. The fact that the optimal ef- 
ficiency of (C/(l -I- C)) 2 represented by the dashed lines 
in Fig.[5]is still achieved by the dotted curves, despite the 
use of truncated controls, proves that truncation of the 
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FIG. 2: Breakdown of optimal adiabatic storage in a cavity 
at TG7 < 10. In (a), the total efficiency of storage followed 
by retrieval is plotted for A = and C = 1, 10, 100, and 1000. 
The horizontal dashed lines are the maximal values (C/(l + 
C)) 2 . Dotted lines are obtained for the input from Eq. (|34[) 
using adiabatic Eq. (|26[l to shape the storage control but then 
using exact Eqs. (I10|l and to numerically compute the 
efficiency. In (b), the same plot is made for C = 10 and 
A/7 = 0,1,10, 100, and 1000. 



storage control does not significantly affect the storage 
efficiency. Since the retrieval field generation is directly 
related to optimal storage by time reversal, as explained 
in Sec. IIV[ this also means that truncating retrieval con- 
trols does not significantly affect the precision with which 
a given retrieval mode e(i) can be generated. The losses 
associated with truncation are insignificant only if the 
conditions in Eq. (fl"8]) and Eq. (|B1[) are still satisfied 
for the truncated retrieval and storage control fields, re- 
spectively. If the limit on the control pulse energy is so 
tight that these conditions are not satisfied, a separate 
optimization problem, which is beyond the scope of the 
present paper, has to be solved. 



VI. FAST RETRIEVAL AND STORAGE 

We have shown that adiabatic storage allows us to 
store optimally a mode of duration T, having any smooth 
shape and any detuning A, provided that the adiabatic- 
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ity condition TC7 ^ 1 is satisfied. In this section we 
solve Eq. (|12[) analytically in the second important limit, 
the so-called "fast" limit, and show that this limit allows 
for optimal storage of a certain class of input modes of 
duration T~ 1/(6*7). 

The fast limit corresponds to the situation when O is 
very large during a short control pulse (|fi| ^ C7 and 
S> |A|), so that we can neglect all terms in Eq. (fT2")) 
except I lips' and S. This corresponds to keeping only 
terms containing f2 on the right-hand side of Eqs. (fTU]) 
and (|11[) and results in undamped Rabi oscillations be- 
tween optical and spin polarizations P and S. One can 
use this limit to implement a "fast" storage scheme, in 
which the input pulse is resonant (A = 0) and the con- 
trol pulse is a short ir pulse at t = T, as well as fast 
retrieval, in which the control is a 7r-pulse at t = T T . 
Provided the 7r pulse is applied on resonance and at ap- 
proximately constant intensity, the condition for perfect 
7r pulse performance is |0| ^S> C7 (assuming C > 1). 
To fully describe these processes, we furthermore need to 
solve Eq. (fP2)) while the control is off, which can also be 
done analytically. 

During fast retrieval, assuming the tt pulse takes place 
at time t — instead of time t = T T and assuming the 
tt pulse is perfect, the initial S — 1 results after the tt 
pulse in P = i. We then solve for P(t) from Eq. (|TU|) and 
insert the solution into Eq. (J5)) to obtain 

£ out (t) = -V2^Ce'^ 1+c ^. (35) 

Consistent with the general expression in Eq. (fT4| and 
the branching ratio argument in Sec. IIV1 the retrieval 
efficiency is again C/(l + C). 

An alternative explanation for why the fast retrieval 
gives the same retrieval efficiency as the adiabatic re- 
trieval is that, thanks to the adiabatic elimination of 
P, the adiabatic limit effectively describes a two-level 
system. Therefore, Eq. ([33)) is in fact a special case of 
Eq. CE1]) with 

n(t) = ( 1 (l + C)+iA)e- zAt . (36) 

Although, at this Q(t), Eq. ([16)) is not a good approxi- 
mation to the actual output field because, for example, 
condition (|29|) is not satisfied, this illustrates the equiv- 
alence of the two approaches. 

Since the control field in fast retrieval is not ad- 
justable (it is always a perfect tt pulse), fast retrieval 
gives only one possible output mode, that of Eq. ([55)) . 
By time reversal, the time reverse of this mode of du- 
ration T ~ 1/(7(7) is thus the only mode that can be 
optimally stored (with efficiency C/(l + C)) using fast 
storage at this C. 

For completeness and to confirm the time-reversal ar- 
gument, the optimal input mode for fast storage can also 
be calculated directly. For an input mode £i n (t) that 
comes in from t = to t = T, assuming a perfect tt pulse 
at t = T, we find by a method similar to the one used in 



fast retrieval that 



S{T) = VTTc J ( 3? ) 

where 

f(t) = -e^ 1+c ^ t - T ^2(l + C) 1 . (38) 

Similarly to retrieval, Eq. ([3"8")) is a special case of Eq. f!M)) 
with n(t) = (7(1 + C) - zA)exp(zA(T - t)) (the time 
reverse of the right-hand side of Eq. ([36)) ). Since f(t) 
is real and normalized according to f(t) 2 dt = 1, this 
integral is a scalar product similar to Eq. ([25)) discussed 
in Appendix [Bj and the optimal fast storage efficiency of 
C/(l+C) is achieved for a single input mode £j n (t) = f(t) 
(up to an arbitrary overall unimportant phase). This 
optimal & m (t) is precisely the (renormalized) time reverse 
of the output of fast readout in Eq. ([33)) , as expected by 
time reversal. 

Two comments are in order regarding the optimal in- 
put pulse £i n (t) = f(t). First, we would like to note 
that short exponentially varying pulses, as in our opti- 
mal solution £i n (t) — f(t), have been proposed before to 
achieve efficient photon-echo-based storage [57j ■ Second, 
it is worth noting that, although £i n (t) = f(t) gives the 
optimal storage, generating such short exponentially ris- 
ing pulses may in practice be hard for high C. Since the 
efficiency is given by the overlap of & m (t) with f{t) (see 
Eq. (|37p ). fast storage in a cavity is inferior in this re- 
spect to fast storage in free space, because in the latter 
case any input pulse satisfying T7 <C 1 and iff 7 ^> 1 
results in storage efficiency close to unity [25j . 



VII. SUMMARY 

In conclusion, we have treated in detail the storage and 
retrieval of photons in homogeneously broadened A-type 
atomic media enclosed in a running- wave cavity. We have 
shown that, provided that no excitations are left in the 
atoms at the end of the retrieval process, the retrieval 
efficiency is independent of the control and the detun- 
ing and is equal to C/(l + C). We have also derived 
the optimal strategy for storage in the adiabatic and fast 
limits and, therefore, demonstrated that one can store, 
with the optimal efficiency of C/(l + C), any smooth in- 
put mode satisfying TC7 S> 1 and having any detuning 
A and a certain class of resonant input modes satisfying 
TC7 ~ 1. We have also noted that the optimal storage 
control field for a given input mode is the time reverse of 
the control field that accomplishes retrieval into the time 
reverse of this input mode. This fact and the equality of 
maximum storage efficiency and the retrieval efficiency 
are, in fact, the consequence of a general time-reversal 
argument to be presented in detail in paper II. In paper 
II, we will also present the full discussion of photon stor- 
age in homogeneously broadened A-type atomic media in 
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free space, while in paper III, we will consider the effects 
of inhomogeneous broadening on photon storage. 

Finally, it is important to note that, to achieve the op- 
timal efficiencies derived in the present paper, it is neces- 
sary to have rigid temporal synchronization between the 
input pulse and the storage control pulse, which may be- 
come difficult in practice for short input pulses. In fact, 
since there is only one accessible atomic mode in the case 
of homogeneously broadened media enclosed in a cavity 
(unless one varies the angle between the control and the 
input |9(), this temporal synchronization is necessary to 
obtain high efficiencies even if the cooperativity param- 
eter is very large. This problem can, however, be alle- 
viated whenever multiple atomic modes are accessible, 
which is the case for homogeneously broadened media in 
free space considered in paper II and for inhomogeneously 
broadened media considered in paper III. In those cases, 
infinite optical depth allows one to achieve unit efficiency 
without rigid synchronization. However, despite this dis- 
advantage of the cavity setup, we will discuss in paper II 
that the cavity setup is superior to the free-space setup 
in other respects, such as the enhancement of the optical 
depth by the cavity finesse and the avoidance of the un- 
favorable scaling of the error as 1/yN (vs 1/C oc 1/N), 
which sometimes occurs in the free-space model. 
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APPENDIX A: DETAILS OF THE MODEL AND 
THE DERIVATION OF THE EQUATIONS OF 
MOTION 

In Sec. IIII1 we gave a brief introduction to the model 
and presented the equations of motion without deriva- 
tion. In this appendix, the details of the model and 
the derivation of the equations of motion Q-© are pro- 
vided. 

The electric-field vector operator for the cavity field is 
given by 

(*) = ex (Jj^) V2 (ae-Wc + a t e -*-i*/c) , (A1) 



where a) is the mode creation operator, U\ is the mode 
frequency, t\ is the polarization unit vector, eo is the 
permittivity of free space, V is the quantization volume 
for the field, and c is the speed of light. 



The copropagating single-mode classical plane-wave 
control field with frequency io 2 is described by an electric- 
field vector 



E2(z, t) = e 2 £2(t) cos(to 2 (t - z/c)), 



(A2) 



where e 2 is the polarization unit vector, and £ 2 {t) is the 
amplitude. Then, using the dipole and rotating-wave ap- 
proximations, the Hamiltonian is 



H = H + V 
H 



N 



frio^a + 2J {tw se a\ a + frw ge a l ee ) , 

i=l 

N 

V = -^d i -(E 2 (z l ,f) + E 1 (z l )) 

N 



(A3) 
(A4) 

(A5) 

+h.c. 



Here h.c. stands for Hermitian conjugate, fr* — 
is the internal state operator of the ith atom between 
states \fi) and \u), Zi is the position of the ith atom, di 
is the dipole moment vector operator for the ith atom, 
f2(t) = i(e|(di • e 2 )\s)i£ 2 (t)/(2h) (assumed to be equal 
for all i) is the Rabi frequency of the classical field, and 

g = i(e|(dj • e i)lg)»^ 2fc%y ( assumecl to be rea l f°r s™ - 
plicity and equal for all i) is the coupling constant be- 
tween the atoms and the quantized field mode. We note 
that, in order to avoid carrying extra factors of 2 around, 
is defined as half of the traditional definition of the 
Rabi frequency, so that a it pulse, for example, takes 
time 7r/(2fi). 

In the Heisenberg picture, we introduce slowly varying 
collective atomic operators 



' sg 



i 
i 
i 

E* 



-IU)2 {t — Zi I c) 



-iuj\ (t — Zi j c) 



-i(oJi— a>2)(*— Zi/c) 



(A6) 
(A7) 
(A8) 
(A9) 



and a slowly varying cavity mode annihilation operator 

£ - ae iuJlt , (A10) 
which satisfy same-time commutation relations 

[a ^(t), a a pit)} = 6 va a^(t) - S^i3a- av (t), (All) 

'£(*),#(*)] = 1, (A12) 

and yield an effective rotating frame Hamiltonian 

H = hAa ee - (Ml(t)a es + hg£a eg + h.c). (A13) 
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The equations of motion are then given by 



£ 



'gs 



—k£ 

Jes&e> 

-(7 + iA)a ge 
iA)a e! 



~ igZ&eg 

- i£lcr es - 



2k£ j n , 
- ig£^ a 



F„ 



, ?e - ig£a, 



ie +ig£a eg 
gs + ig£ (tgg ' 
)- 



- ig£ ] a ge + F ee , 

ig£^a gs + F es , 
(A14) 



with the input-output relation for the quantum field 
given by Eq. © . 

In Eqs. (|A14[) we have introduced decay, which, in turn, 
necessitated the introduction of Langevin noise operators 
F^ for the atomic operators and the input field £- lVL for 
the quantum field. The radiative decay rate of the ex- 
cited state |e) is ~f c = 7 0S +7cg, the sum of decay rates into 
\s) and into \g). The decay rate of optical coherence a ge 
is 7 = 7o/2 + 7dcph where, in addition to radiative decay, 
we allow for extra dephasing, such as, for example, that 
caused by collisions with buffer gas atoms in warm vapor 
cells. Similarly the decay rate 7' = 7 c /2 + j' dcph of a es 
allows for possible extra dephasing, while the decay rate 
7 S of & gs is due to dephasing only. In some experiments 
@ , 7 S comes from the transverse diffusion of atoms out of 
the region defined by the quantized light mode. In these 
cases, the decay of a gs will be accompanied by popula- 
tion redistribution between states \g) and \s). In order 
to ensure that the corresponding incoming noise is vac- 
uum (which our analysis requires, as we explain below 
and in Sec. |TT|, we will assume in such cases that the in- 
coming atoms are fully pumped into the level \g), which 
would correspond to a 27 s decay rate of a ss into a gg (not 
included in Eqs. (I A14[l since it does not affect the final 
equations). This is indeed the case if, as in Ref. @, the 
control beam diameter is much greater than the diameter 
of the quantized light mode. 

Assuming that almost all atoms are in the ground state 
at all times (a gg 
polarization P = a 
and working to first order in £, we obtain Eqs. p ]) -([S ]) . 
where F P = F ge /yfiqN and F s = F gs /y/2^N. 

Using the generalized Einstein relations (EH, [H^] 

(F^(t)F al3 (t')) = {D(a^a a p) - D^^&ap 

-a^Di&^Sit-t'), (A15) 

where D(a^) denotes the deterministic part (i.e., with 
noise omitted) of the equation for in Eqs. (|A14|) . and 
again using the approximation that almost all atoms are 
in the ground state, we find that the only nonzero noise 

>t 
p- 

(Fs(t)Fl(t')) 



N and a ss w a ee ps a es ps 0), defining 
ge/V^V and spin wave S = a gs /\/~N, 



correlations between Fp, Fs, Ft, and Ft an 



(F P (t)F P (t')) 



8{t-t'). (A16) 



The fact that normally ordered correlations are zero 
means that the incoming noise is vacuum, which is pre- 
cisely the reason why, as noted in Sec. |IIJ efficiency is 



the only number we need in order to fully characterize 
the mapping. The property of Eqs. (|A14p that guaran- 
tees that the incoming noise is vacuum is the absence of 
decay out of state \g) into states |e) and \s). The decay 
into state |e) does not happen because the energy of an 
optical transition (on the order of 10 4 K) is much greater 
than the temperature, at which typical experiments are 
done. In contrast, the energy of the \s) — \g) transition 
in some experiments, such as the one in Ref. 0], may 
be smaller than the temperature. However, the \s) — \g) 
transition is typically not dipolc allowed, so that the rate 
of \g) decay into \s) can be neglected, as well. As noted 
above, for the case when atoms are diffusing in and out 
of the quantized light mode, to keep the decay rate of 
\g) zero, we assume that the incoming atoms are fully 
pumped into |<7). 

From Eq. (|A11[) and with the usual a gg ~ N assump- 
tion, we have 



>(t),pt( t ) 



= 1, 



(A17) 
(A18) 



In particular, this means that S can be thought of as an 
annihilation operator for the spin-wave mode, into which 
we would like to store the state of the incoming photon 
mode. 

The input and output fiel ds, which propagate freely 
outside of the cavity, satisfy [5J] 



&ut(t),£L(f) =5(t-t') (A19) 



and can be expanded in terms of any orthonormal set 
of field (envelope) modes {h a (t)} defined for t € [0, 00), 
satisfying the orthonormality relation L dth* a {t)h^{t) = 
6 a p and completeness relation ^ Q h* a (t)h a (t') = S(t—t'), 
as 



£ia(t) = 2J h a (t)a a , 
a 

£out(*) = X! h a (t)b a , 



(A20) 
(A21) 



where annihilation operators {a a } and |^q| f° r the in- 
put and the output photon modes, respectively, satisfy 



(A22) 



Repeating for clarity the setup from Sec. IIIH we recall 
that all atoms are initially pumped into the ground state, 
i.e., no P or S excitations are present in the atoms. We 
also assume that the only input field excitations initially 
present are in the quantum field mode with annihilation 
operator clq and envelope shape ho(t) nonzero on [0, T]. 
The goal is to store the state of this mode into S and 
at a time T T > T retrieve it back onto a field mode. 
During storage, we can, in principle, solve the operator 
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Eqs. I©-©) for S(T) as some linear functional of £;„(i), 
F P (t), F s {t), 5(0), and P(0). The storage efficiency is 
then given by 

(number of stored excitations) (S' t (r)S'(r)) 



(number of incoming photons) £dt{E} n (t)£ in (t)) ' 

(A23) 

Since S(0) and P(0) give zero when acting on the initial 
state, and since all normally ordered noise correlations 
are zero, only the term in S(T) containing £ m (t) will 
contribute to the efficiency. Moreover, ho(t)ao is the only 
part of £i n (i) that does not give zero when acting on the 
initial state. Thus, for the purposes of finding the stor- 
age efficiency, we can ignore F P and F s in Eqs. ©-© 
and treat these equations as complex number equations 
with P(0) = 0, S{0) = 0, and £ ln (t) = h Q (t). We have 
here dropped the carets on the operators to denote their 
complex number representations. To get back the non- 
vacuum part of the original operator from its complex 
number counterpart, we should just multiply the com- 
plex number version by clq- 

Similarly, during retrieval, we can ignore Fp(t) and 
Fs(t) and can treat Eqs. I©-© as complex number equa- 
tions with the initial and boundary conditions given in 

secnni 



APPENDIX B: SHAPING THE CONTROL FIELD 
FOR THE OPTIMAL ADIABATIC STORAGE 

In this appendix, we present the derivation of Eq. (|2U1) . 
which gives the optimal storage control field during adi- 
abatic storage. We then verify that this optimal control 
is just the time reverse of the control that retrieves into 
the time reverse of the input mode. 

To solve for the control field O(i) that maximizes the 
storage efficiency t] s in Eq. (j2l)|) , we note that f(t) defined 

in Eq. satisfies f Q T \f(t)\ 2 dt < 1, with the equality 
achieved when 



2 7 (l + g) 
7 2 (1 + C) 2 + A 



rM0,T)»E 



(Bl) 



which is equivalent to the requirement we had in Eq. fT8J) 

for complete retrieval. Since we also have L \£i n (t)\ 2 — 
1, the integral in Eq. (|25|) can be seen as a simple scalar 
product between states, and the efficiency is therefore 
Vs < C/(l + C) with the equality achieved when (up to 



an undefined overall unimportant phase) 



(B2) 



We will now show that, for any given £i n (t), A, and 
C, there is a unique control that satisfies Eq. (|B2[) and 
thus give s the maximum storage efficiency C /(l + C). In 
Refs. [60l.l6l|. this control was found through a quantum 
impedance matching Bernoulli equation obtained by dif- 
ferentiating Eq. (|B2[) . In order to be able in paper II to 
generalize more easily to free space, we will solve Eq. (|B2[) 
directly. To do this, we follow a procedure very similar 
to that in Sec. IV Al We integrate the norm squared of 
Eq. (HH from to t to get 



2h(t,T)-y(l + C) 



2h(0,T) T (l + C) 



\& m (t')\ 2 dt' = e ^(i+cja+A* _ e ^T+W+Z* . (B3) 



Since h(T, T) = 0, the normalization of & m (t) requires the 
satisfaction of Eq. (jBlj) . Assuming it is satisfied to the 
desired precision, we solve Eq. (|B3p for h(t,T), and then 
taking the square root of the negative of its derivative 
with respect to t, we find |0(t)|. Knowing h(t,T), the 
phase of f2(t) can then be determined from Eq. (|B2[) . 
Putting the magnitude and the phase together, we obtain 
the expression for the optimal control given in Eq. (|26p . 

We will now show that, as expected from the time- 
reversal argument, the optimal control we derived is just 
the time reverse (0,(t) — > 51* (T — t)) of the control that 
retrieves into £* n {T — t), the time reverse of the input 
mode. To see this, we note that from Eq. (f2T)]) it follows 
that the magnitude of the control field f2(t) that retrieves 
into e(t) = £* n (T — t) is determined by 

dt'\£ in {t')\ 2 = e -^W+^ /l(0 ' T - t) . (B4) 

Putting e(t) = £* n (T - t) into Eq. (|2T|). taking the com- 
plex conjugate of the result, and evaluating at T — t, we 
get 



Q*(T-t)-- 



7 (l + C)-iA 



£in{t) 



■ Afo(0,T— t) 



dt' 



(B5) 

Since h(0,T - t) determined by Eq. jB4| and h(t,T) 
determined by Eq. (|B3[) are equal, the right-hand side 
of Eq. (|B5j) is, in fact, equal to the right-hand side of 
Eq. (|2l)|) , as desired. 
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